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^ ■ 1 INTRODUCTION 



ABSTRACT 

We present a halo model prediction of the image separation distribution of strong 
lenses. Our model takes account of the subhalo population, which has been ignored 
in previous studies, as well as the conventional halo population. Halos and subhalos 
are linked to central and satellite galaxies by adopting an universal scaling relation 
between masses of (sub-)halos and luminosities of galaxies. Our model predicts that 
10%— 20% of lenses should be caused by the subhalo population. The fraction of lensing 
by satellite galaxies (subhalos) peaks at ~ 1" and decreases rapidly with increasing 
image separations. We compute fractions of lenses which lie in groups and clusters, 
and find them to be ~ 14% and ~ 4%, respectively: Nearly half of such lenses are 
expected to be produced by satellite galaxies, rather than central parts of halos. We 
also study mass distributions of lensing halos and find that even at image separations 
of ~ 3" the deviation of lens mass distributions from isothermal profiles is large: At or 
beyond ~ 3" image separations are enhanced significantly by surrounding halos. Our 
model prediction agrees reasonably well with observed image separation distributions 
from galaxy to cluster scales. 

Key words: cosmology: theory — dark matter — galaxies: elliptical and lenticular, 
cD — galaxies: formation — galaxies: halos — galaxies: clusters: general — gravita- 
tional lensing 



X 



Image separations between multiple images play an impor- 
tant role in the statistics of lensed quasars: The image sep- 
aration is mainly determined by the potential wall depth 
of the lens object, thus the image separation distribution 
reflects hierarchical structure in the universe. This fact sug- 
gests that lensed quasars may categorized roughly into two 
populations; small-separation lens (8 ~ 1") that are pro- 
duced by a single galaxy, and a large-separation lens [9 > 
10") that are caused by clusters of galaxies. Thus far about 
80 lensed quasars are known: Most are sma ll-separation 
lense s and only one len s, SDSS J1004+4112 lllnada et alJ 
120031: lOeuri et alJ l2004h . has the image separation larger 
than 10". The current largest homogene ous lens survey, the 
Cosmic Lens All S ky Survey (CLASS: iMvers et aLlEooll 
iBrowne et al"ll2003h . contains small separation only, thus it 
may not be suitable for studying the image separation dis- 
tribution from galaxy to cluster scales. However, ongoing 
lens surveys such as that in the Sloan Digital Sky Survey 
(SDSS) are expected to unveil the full distribution, which 
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will be extremely useful for understanding the assembly of 
structures. 

In theoretical sides, there have been many attempts to 
calculate the full image separation distribution. Historically, 
the image separation distributio n is computed from either 
the galaxy luminosity func tion (farmer et al.l | l984l) or the 
mass function of dark halos iNaravan fc Whitejl988l) . While 
the former approach accounts for the observed lensing prob- 
ability and its distribution at 9 < 3" fairly well, it has diffi- 
culty in explaining the existence of the large-separation lens 
for which dark matter dominates the lens potential. There- 
fore in this paper we take the latter approach. 

However, modeling the image separation distribution on 
the basis of dark halos is not simple. It has been found that 
none of models that consider only one population for lens- 
ing halos are consistent with current observations. Correct 
models need to include a characteristic mass around Af coo i ~ 
10 13 Mq , where the density distribution inside dark halos 
changes JPorciani fc Madaull200ol: iKochanek fc Whitel200l| ; 
KeetorJl200ll: ISarbu et alJl200li iLi fc Ostrikerll2002l booa 
Qguri l20oJjMal |2003t iHuterer fc Mal|2004 iKuhlen et all 
2004 IChenl 12004 : IZhand 12004 IChenl 120051 : IZhang et all 
2005): Large mass halos M > M coo \ (correspond to groups 
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or clusters) have cooling time longer than the Hubble time, 
thus the internal structure of dark halos is not strongly af- 
fected by bar yon cooling and is w ell described by that of 
dark matter iNavarro et al.l [l997fl . On the other hand, in 
small mass halos M < M coo i (correspond to galaxies) the 
baryon cooling is very efficient and the mass distribution be- 
comes strongly centrally concentrated, which is well approx- 
imated by the Singular Isothermal Sphere (SIS). Although 
this two-component model is successful in explaining the 
image separation distribution to some extent, clearly it is 
not sufficient. First of all, this model considers only galax- 
ies that lie at the center of isolated dark halos, and thus 
does not take galaxies in groups and clusters into account. 
Put another way, this model neglects substructures (sub- 
halos) which should corresponds to satellite galaxies. Actu- 
ally, it has been speculated that lens systems in suc h dense 
environments is quite common jKeeton et al]|2000t) . In ad- 
dition, sometimes lens systems are significantly affected by 
the mass distributio n outside the lens object. A good ex- 
ample is Q0957+561 JWalsh et al.lll979l) : its lens object is a 
galaxy in a cluster, and the image separation is much larger 
than expected from the luminosity of the lens galaxy because 
the image separation is boosted by the cluster potential. In- 
deed, it has been shown theoretically that at intermediate- 
separation regime (8 — 3" — 7 ") lens galaxies in dense envi- 
ronments are quite common dOguri et alJl2005l) . Therefore 
the correct model needs to account for such external mass. 

To construct a model that is based on dark halos and 
subhalos, we have to model the relation between mass of 
dark halos (subhalos) and luminosities of galaxies inside 
them. There are several models to link them: Conditional 
luminosity function (CLF) which is defi ned by the lumi- 
nosity function in a halo of given mass llYang et alJl2003t 
van den Bosch et al]l2005T) . universal mass- luminosity rela- 
tion b etween halo/subhalo m ass and hosted galaxy lumi- 
nosity llVale fe Ostrikeill2004t) . and the mixture of these tw o 
models dCoorav fe Milosavlievicteo'oll ICoorav fe Cenl2005l) . 
There models are calibrated to match observed galaxy lu- 
minosity functions, clustering of galaxies as a function of 
galaxy luminosities, and the luminosity function of galaxies 
in groups or clusters. 

In this paper, we study the image separation distribu- 
tion with a particular emphasis on the different contribu- 
tions from central and satellite galaxies. For central galaxies, 
we consider baryon cooling i n a d ark mat t er ha lo as done 
bv iKochanek fe White! (1200 ill and iKeetonl (1200 ll) . Satellite 
galaxies, which have been neglected in previous studies, are 
linked to dark halo s ubstructures. We basica lly take the ap- 
proach proposed bv lVale fc Ostrikerl J2004) to relate mass 
of dark halos and subhalos with luminosities of galaxies, but 
we modify it to account for the relativ ely small abundance o f 
satellite galaxies in galactic halos fsee lCoorav fc Ccn 2005). 
For satellite galaxies, the external mass from the host halo, 
which is shown to be important at large-image separations, 
is taken into account. In addition to the difference between 
central and satellite galaxies, we also pay particular atten- 
tion to how lens objects make the transition from SIS to 
NFW as halo masses increase. 

This paper is organized as follows. In SJ5J we introduce 
models of dark halos and subhalos that are used to compute 
lensing probabilities. Section [3] are devoted to model the re- 
lation between (sub-)halo masses and galaxy luminosities. 



We compute the image separation distributions in S0] and 
show the results in SjS] Finally, we summarize and discuss 
our results in [JS] Throughout the paper we assume a con- 
cordance cosmology with the mass density Om = 0.3, the 
cosmological constant = 0.7, the dimensionless Hubble 
parameter h = 0.7, and the normalization of the mass fluc- 
tuation erg = 0.9. 



MODELING DARK HALOS AND 
SUBHALOS 



2.1 Dark Halos 



For the density profile of original dark halo s, we adopt the 
spherical NFW profile (iNavarro et alJll997l) : 



p{r) = 



(r/r s )(l + r/r 



(1) 



The characteristic density p B is computed from the nonlin- 
ear overdensity A v i r - The scale radius r s is related to the 
concentration parameter c = r v i r /r s , where r v i r is the virial 
radius. The median concentration parameter depends on the 
mass M and red shift z of the dark ha lo: We adopt the fitting 
form derived bv lBullock et alJ J200ll) : 
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where M* is determined by solving om — S c at z = (om 
is the standard deviation of linear density field smoothed 
with a top-hat filter enclosing mass M, and 5 C ~ 1.68 is 
the threshold linear overdensity above which the region col- 
lapses) . We also include the scatter in the concentration pa- 
rameter, which is well described by a log-normal distribu- 
tion: 



p(c) = 



2-Ka h, 



■ exp 



(lnc- In cf 



2a, 2 

In c 



(3) 



where the standard deviation is <ri nc = 0.3. 

In this paper, we consider the modification of mass 
distribution in dark halos due to baryon cooling. Tra- 
ditionally, the response of dark matter to baryon infall 
has been calculated by the model of adiabatic contraction 
(Blu menthal et al][T9 86.). The model assumes the spherical 
symmetry and computes the response by imposing the con- 
servation of angular m omentum between b efore and after 
baryon infall. Recently, iGnedin et al.l (|2004) studies the va- 
lidity of the model using high-resolution numerical simula- 
tions, and found that the model can be improved by taking 
the eccentric orbits of particles into account. They also pre- 
sented a series of analytic fitting functions which make it 
much easier to implement the modification. We adopt this 
modified adiabatic compression model to derive the total 
mass distribution. For the final mass d istribution of the 
cooled baryonic component, we assume iHernauisJ (Il99d) 
profile: 



Pb{r) = 



2tt (r/r b )(r b + r) 3 ' 



(4) 



which is known to describe the stellar density profile of ellip- 
tical galaxies well. Note that to compute the final mass dis- 
tribution we need to specify the cooled mass Mb (or equiva- 
lently mass fraction /b = Mb /M) and the core radius rb (or 



© 2004 RAS, MNRAS OOO.ITlliTI 



The image separation distribution of strong lenses 3 



equivalently the ratio of the core radius to the scale radius 
Xb = Tb/r s ), in addition to halo parameters such as c and 
M. 

For the mass function of dark halos, we use the form 
proposed bv lSheth fc Tormenl il999T> : 
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■ exp 
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(5) 
(6) 



where p(z) is the mean matter density at redshift z, A = 
0.322, a = 0.7 07, and p = 0.3. The m ass function is the 
modification of lPress fc Schechted Jl974T) to account for el- 
lipsoidal nature of gravitational collapse in cosmological sit- 
uations, and is known to reproduce the mass distribution 
and its redshift evolution in numerical simulations (e.g., 
iReed et aljl2003Tl . 



2.2 Subhalos 

We consider subhalos which lie in a dark halo with mass M. 
The mass of the subhalo is denoted by m s and the distance 
from the halo center by r. We always adopt an SIS for the 
density distribution of galaxies associated with subhalos for 
simplicity. Indeed it has been known that the central part of 
the density profile of total (dark matter plus cooled baryon) 
matter in a dark halo after baryon infall is quite close to 
an SIS. For the mass and spati al distributions of s ubhalos, 
we use a model constructed bv lOeuri fc Led i2004T) . In the 
paper, they showed that the analytic model reproduces well 
the results of high-resolution numerical simulations. In real- 
ity, we use the following fitting formulae of mass and spatial 
distributions to speed up the computation. 

The fitting functio n of the number distrib ution of sub- 
structures was given bv lVale fc Ostrikerl J2004t) : 



dN 



0.18 



( 1713 \ 

dm s -yT{2-(3) Um) 



exp 



m 3 \ 1 
"7M J 7M' 



(7) 



where (3 = 1.91 and 7 = 0.39. The mass function has a 
power-law form at m s <C M, and has a sharp cut-off at 
m s ~ M. In addition, we assume that dN/dm s = at 
m s > M. We compare t his fitting function with the re- 
sult of lQguri fc Led <l2004l) in Figure [5] confirming the good 
accuracy of the fitting function. 

It has been shown that the radial number density distri- 
butions of subhalos are less centrally concentrated relative 
to the dark matter particle components. Moreover, they de- 
pend on masses of subhalos: More massive subhalos are pref- 
erentially located in the external regions of their host halo s 
<De Lucia et al.ll2004t lOeuri fc Ledl2004 iGao et allEiol . 
We find that the spatial distribution of subhalos is described 
by an NFW profile (eq. 0) with different concentration pa- 
rameter (which depend on m 3 ) from that of the host halo. 
Specifically, the concentration parameter of the radial dis- 
tribution of subhalos is fitted by 



1 + 



( 3m s /M 
\ 10- 7 



(«) 



where c is the concentration parameter of the correspond- 
ing host halo. Figure [5] shows the comparison of this fitting 
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Figure 1. The comparison of the cumulative nu mber distribu- 
tions of subhalos derived from an analytic model of lQguri fc Leel 
12004) with the fitting function (obtained by integrating equation 
171 over m s ). The former is shown by histograms, and the latter 
is plotted by a solid line. We plot the distributions of the ana- 
lytic model for three different masses of the host halos. Note that 
the cumulative mass distribution derived from the fitting function 
does not depend on the mass of the host halo when it is plotted 
as a function of m s /M. 



with the model of lOeuri fc Led lj2004l) . They agree reason- 
ably well with each other for a range of parameters we are 
interested in. 



3 LINKING DARK HALOS AND SUBHALOS 
TO GALAXIES 

We consider a central galaxy with luminosity L which resides 
in a dark halo with mass M. We adopt the mass-lumin osity 
(bj band) relation proposed bv lVale fc Ostrikeil i2004l) : 



L(M) = 5.7 x 10 9 /i" 2 L Q 



q + M<r r)s 



1/s 



(9) 



where p = 4, q = 0.57, r = 0.28, s = 0.23, and Mu = 
M I (10 11 h~ x Mq) . Basically this relation asymptote to I a 
M 4 at low mas s end and L oc A/ ' 28 at high mass end. 
IVale fc Ostrikerl <|2004|) claimed that the same relation can 
be applied to dark halo substructures (with mass m s ) and 
satellite galaxies (with luminosity L s ): 



5.7 x 10 9 /i" 2 L Q 



(p — r)i 

q + 



l/s 



(10) 



where mu = m/(10 11 h^ 1 Mq) and m is the mass of subhalos 
before the tidal stripping of the outer part of subhalos. We 
adopt the a pproximation relation m = 3m s , which was also 
assumed in IVale fc Ostrikerl <|2004F) . throughout the paper 
for simplicity. 
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Figure 2. The cumulative radial distributions of subhalos in an analytic model iQguri fc Leel2004l) are compared with those of an NFW 
fit with concentration parameters given by equation JgJ. The former distributions are plotted with histograms, while the latter are shown 
by curves. From left to right panels, we change the mass of the host halo. Different lines indicate different masses of subhalos. 



While IVale fc Qstrikel (2004) adopted monotonic one- 
to-one correspondence bet ween halo/subhalo mass e s and 
resident galaxy luminosities, ICoorav fc Milosavlievid (I2005T) 
showed that it is essential to assign log-normal scatter to 
the mass-luminosity relations to recover the tail of the ob- 
served luminosity function. Following them, we assume the 
log-normal luminosity distribution of 



p(L\M) 



27rln(10)EL 



exp 



log 10 [L/L(M)f 
2E 2 



(11) 



with E = 0.25, and the same distribution with mean of 
L s (m s ) for L s . 

Although this model is successful in explaining many 
observations such as the occupation number, the luminos- 
ity function of clusters and the group luminosity function, 
it cannot reproduce the strong dependence of the faint 
end slope of the CLF on the halo mass: Observationally, 
it changes from ~ —1.2 at the cluster scale to ~ at 
galactic scales. Since the mass function of subhalos at low 
mass is dn/dm 3 tx m' 13 where /3 ~ 1.9 (see H2.2L the 
mass-luminosity relation L s (m a ) tx ml, yields the CLF 
oc LJ 1-0 ' 9 ^ . Therefore the above mass-luminosity relation 
implies the faint end slope of ~ —1.2. It is also found that 
even if we modify the value of r\ we cannot make the faint end 
slope shallower than —1. Hence to account f or the observed 
trend we need extra mechanism. Following ICoorav fc Cenl 
(2005), we introduce an efficient function to suppress the 
number of satellite galaxies in galactic halos: 



f(m a ) = 0.5 



1 + erf 



log m — log m c 



(12) 



where <j m = 1.5 and the cutoff mass m c depends on the host 
halo mass M : 



m c = 10 11 h 1 M(? ) exp 



AI 



5 x lO^/i-iM, 



(13) 



For less massive (M < 10 13 fc _1 M Q ) halos, the efficiency 
function behaves as f(m a ) — > when m s <C 10 n /i -1 M G 



and f(r 



1 when m s 3> 10 11 h 1 Mq. At cluster scales 



M > 10 14 h" 1 Mq the efficiency function is almost always 
unity. 

We neglect the redshift evolution of the relations pre- 
sented in this section, since most lens galaxies are at z < 1 
where the evolution is expected to be mild. 



4 COMPUTING THE IMAGE SEPARATION 
DISTRIBUTION 

To compute lensing cross sections, we have to convert galaxy 
luminosities to any quantities which characterize the mass 
distribution of galaxies. To do so, we use the scaling relations 
between the galaxy luminosity L, velocity dispersion <r, and 
effective radius i?o- Specifically, we ad opt those derived from 
9,000 early-type galaxies in the SDSS dBernardi et all2003T) : 



l T0 Liq 



10 102 /if 2 L Q V 10 °- 52/l 7o lk P c 



Bo 



(14) 
(15) 



where L is (/-band galaxy luminosity and we defined ft.70 = 
h/0.7. Note that the effective radius is related to the core 
radius of Hernquist profile (eq. [1]) by rb = 0.551i?o- For 
simplicity, we neglect tiny difference of 6j-band and g-band 
luminosities and regard L in the above equations as the same 
luminosities as those in Although the fundamental plane 
can describe the correlation more tightly, we adopt these 
relations between two observables because they are much 
easier to handle. We also convert luminosities of galaxies to 
stellar masses assuming the universal stellar mass-to-light 
ratio of T = 3.O/17OM0/LQ. 

One of the most important ingredient in computing 
lensing probabilities is magnification bias which is deter- 
mined by the source luminosity function. Below we assume 
a power-law <3?(L) oc , which makes the computation of 
magnification bias easier. 
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4.1 Dark Halo Component 

As discussed in H2.ll we use the im proved adiabatic compres- 
sion model of lOnedin et al.l 1)200 4) to compute the total mass 
distribution of a galaxy plus dark halo and hence to derive 
its lensing cross section. First, we derive the luminosity of a 
given halo with mass M using the mass-to-luminosity rela- 
tion ©, and estimate the cooled baryon fraction /b and the 
normalized core radius Xh from the stellar mass-to-light ratio 
T and equation 1151 . respectively. From the mass distribu- 
tion, we compute an image separation 6 and bias ed cross sec- 
tion g iens: We use approximations proposed bv lOeuri et alJ 
(2002) to compute these quantities. Then the image sepa- 
ration distribution can be obtained by integrating the cross 
section: 



dP h f, n , x3«2i 



dM — a lcnE 5(6-6(M)). (16) 



The final probability distribution is obtained by averag- 
ing equation 11611 over the probability distribution functions 
(PDFs) of the concentration parameter and the mass-to- 
luminosity relation. 

4.2 Subhalo Component 

The projected mass distribution of a subhalo is approxi- 
mated as the sum of an SIS and external convergence K ox t 
originates from the mass associated with the host halo. We 
only consider the external convergence and neglect the ef- 
fect of the external shear, since it has been shown that the 
external convergence is mor e effective to mod ify the im- 
age separation distribution JOeuri et al.ll2005l) . From the 
velocity dispersion a, the image separation 8 and the bi- 
ased cross section cri en s c an be computed analytically (see 
iKeeton fc Zabludofl1l2004l . for the dependence on K ox t): 



20E 



(1 — k c 



(Tie 



7T(D 



87r(q/c) 2 As 
(1 — K cxt ) D os ' 

a 2*(1 - /w)- 2 ^ 
3-£ 



(17) 



(18) 



where D os is the angular diameter distance from the observer 
to the source, etc. The external convergence is computed 
from the mass distribution of the host halo after baryon 
cooling f H2.HI . Since the external convergence is written as 
a function of the projected radius R, the PDF of the external 
convergence becomes 



p(/Sext)d/s ex t = fn(R) ^j™* dR, 
dR 



(19) 



where fn(R) is the PDF of subhalos projected along line- 
of-sight, which we adopt the projected NFW profile with its 
concentration parameter given by equation (S): 



fn(R)dR 



r/2 
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( r vil ) 



'9 I — r L — dR, 

smro / r v ir 



(1 + z) 2 
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(21) 
(22) 



h{x) 



dyg(y) 
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i 



-1 + 



V l—x l - 



= arctanh 



l-x 
1 + x 



= arctan 



(x > 1) 
(x > 1). 



(23) 



From these, we compute the image separation distribution 

as 

dPs f , ,., , 3 cdt f ,,,dn f , . dN 

W = dZl{1+Zl) M dM dM dmsf{ms) d^s 



X / C^K ex tp(K ex t)o■lens'5(6 , - (/Start))- 



(24) 



As in H4.ll the final probability distribution is obtained by 
averaging equation 1241 over the PDFs of the concentration 
parameter of host halos (eq. [3]) and the mass-to-luminosity 
relation (eq. 1111 ). We set the upper limit of the integral over 
«: cx t to 0.9, because in the events beyond the limit (« e xt ^ 

0. 9) subhalos should fall inside the critical radius of the host 
halo and therefore they should be included in lensing by dark 
halos computed in H4.ll 

We note that the inclusion of the efficient function (eq. 
1121 ) implies that a part of low-mass subhalos remains dark, 

1. e., does not harbor a galaxy. We neglect lensing by these 
dark subhalos because their contribution to total lensing 
probability distributions is ra ther small unless the inner pro - 
file of subhalos is very steep (lMall2003l: iLi fc Ostrikeill2003D . 

The total image separation distribution is simply given 
by a sum of equations 11611 and 1241 : 

dP t _ dPj, dP s 

dO d6 dO ' ( ' 



5 THE IMAGE SEPARATION DISTRIBUTION 

In this section, we compute image separation distributions 
for halos (central galaxies) and subhalos (satellite galaxies). 
In the specific examples below, we fix the source redshift to 
z s = 2.0 which is typical for lensed quasars. The magnifi- 
cation bias is computed assuming a power-law luminosity 
function with £ = 2.1 which is the lu minosity function of 
the CLASS survey llMvers et al] 120031) . First we will show 
the image separation distributions for dark halos and sub- 
halos separately, and next we see the contribution of each 
component to the total distribution. 



5.1 Dark Halo Component 

Before going to the image separation distribution, we check 
the image separation as a function of halo masses to see 
the effect of baryon infall on lensing, which is shown in Fig- 
ure |S| Two special cases are also shown for reference: One 
is neglecting baryon infall and leaving halos as NFW, and 
the other is adopting an SIS approximation with the ve- 
locity dispersion inferred from the luminosity of the central 
galaxy. It is found that the image separations agree well with 
those of SIS at low-mass ends, which is consistent with ob - 
servations fe.g.. lRusin et al.l2003tlRusin fc Kochanekl2005l) . 
With increasing masses, image separations begin to deviate 
from SIS from 6 ~ 1" (M ~ to 12 /! -1 M ). At 6 > 10" 
(M > 3 x 10 14 /i _1 A/q) image separations are quite simi- 
lar to those of NFW. This figure justifies to some extent 
the two-population model in which SIS and NFW lenses are 
considered, but we note that at M ~ 10 13 - 10 14 /i _1 M o 
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Figure 3. The image separations of lensing by dark halos after 
baryon cooling are plotted as a function of halo masses. In this 
plot, we do not take the PDFs of the mass-to-luminosity relation 
and concentration parameter into account. The lens redshift is 
fixed to zi = 0.5. For comparison, we also show cases that we 
neglect baryon infall (NFW; dashed line) and we adopt an SIS 
approximation for the mass distribution of the halos after baryon 
infall (SIS; dotted line). Note that the velocity dispersion for SIS 
is derived from the scaling relation of equation 1141 . 
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Figure 4. The image separation distribution from the dark halo 
component (eq. 1161 ). Lines are same as Figure 131 



things are too complicated to be described by such simple 
two-population model: Central galaxies are mainly responsi- 
ble for lensing, but halo masses outside the lensing galaxies 
boost the image separation significantly, making the SIS ap- 
proximation inaccurate. 

We show the image separation distribution in Figure [I] 
As in Figure [3] at small- and large-separation regions the 
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Figure 5. The probability distribution of "effective" external 
convergence K a g (eq. 1261 ) which is a measure of the deviation 
of the halo mass distribution from an SIS approximation of the 
corresponding central galaxy. We derive the distributions for sev- 
eral representative image separations: 9 = 0'.'3 (upper left), 1" 
(upper right), 3" (lower left), and 6" (lower right). 



probabilities are close to those of SIS and NFW, respec- 
tively, but at 3" < 9 < 10" it is difficult to describe lensing 
probabilities as either SIS or NFW alone. The lensing prob- 
abilities also differ from the sum of those of SIS and NFW. 
Therefore, we conclude that the image separation distribu- 
tion cannot be described by a simple two-population model: 
We need to model the transition of halo mass distributions 
from SIS to NFW to predict correctly the distribution from 
small- to large-image separations. 

We explore the transition from SIS to NFW more in 
terms of an "effective" external convergence defined by 



w = 1 — 



#SIS 



(26) 



where for each halo #sis is computed from the luminosity 
of the central galaxy via equation l|14|l . The case K e g — 
means that the Einstein radius of the halo is exactly the 
same as that inferred from the luminosity of the central 
galaxy. The definition is motivated by the description of the 
image separation in the SIS plus external shear model (eq. 
1171 ). We compute the PDFs of /t e ff for several represen- 
tative image separations, which are shown in Figure^] As 
expected, at 6 < 1" lens halos are on aver age close to an SIS , 
being consistent with ob servations (e.g.. lRusin et ai1l2003l : 
iRusin fc Kochanekl feoOS). However, at 9 > 3" they clearly 
have larger image separations than those of corresponding 
SIS lenses. Therefore this Figure represents another case for 
the deviation from simple SIS lenses at that image separa- 
tion region: At or beyond 9 > 3" the image separations are 
enhanced significantly by surrounding dark halos. 

We also point out that K cx t has rather broad distri- 
butions even at small-separations 9 < 1", which comes 
from the scatters of the concentration parameter (eq. [3]) 
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Figure 6. The image separation distribution from the subhalo 
component (eq. 1241 ) is plotted by a solid line. We also calcu- 
late the distribution without including re ex t , which is shown by a 
dotted line. 



and the mass-to-luminosity relation (eq. 1111 ) around their 
medians. This scatter around SIS approximations may be 
able to explain possible heterogeneous nature of lens galax- 
ies inferred from the combination of gravitati onal lensing 
and stellar dynamics llTreu fc Koopmandl2004|) or time de- 
lay measurements betw een multiple images of lensed quasars 
llKochanek et al1l2005t) . 



5.2 Subhalo Component 

We show the image separation distribution from the subhalo 
component in Figure |S| We find that the distribution has a 
peak at around 8 ~ 0.5", and rapidly decreases with increas- 
ing 9. It is also seen that the shape of the curve changes at 
around 9 ~ 5" over which the probability decreases less 
rapidly. We attribute this behavior to the contribution of 
very large external convergence («; cx t > 0.5), because the 
probabilities are dominated by such extreme events (see be- 
low). This means that the result is somewhat sensitive to 
the upper limit of /t ext which we set 0.9. However, as we 
will see in the next subsection, this does not cause any sig- 
nificant problems since at that image separation region the 
fraction of subhalo lensing in the total lensing probabilities 
are rather minor. 

As shown in lOeuri et alJ (l2005f) , the external conver- 
gence affects the image separation distribution significantly, 
particularly at the tail of the distribution 9 > 3". To show 
this explicitly, in Figure|5]we also plot the distribution with- 
out including /text- It is found that the external convergence 
enhances the lensing probabilities at 9 > 1" where the prob- 
abilities are de creasing func t ion of 9. This is consistent with 
the finding bv lOeuri et alJ ((2005) . The enhancements are 
quite significant; factors of ~ 6 at 9 = 2" and ~ 80 at 
9 = 3". Note that the s e num bers are much larger than the 
results of lOeuri et alJ (1200 5). because here we concentrate 
on satellite galaxies that should have larger external conver- 
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Figure 7. The distribution of external convergence re ex t f° r the 
subhalo population. The PDFs for several different image separa- 
tions are shown: 9 = 0'.'3 {upper left), 1" {upper right), 3" {lower 
left), and 6" {lower right). 



gence on average than central galaxies. Our result indicates 
that it is very important to include the external convergence 
of satellite galaxies, which originates from the host halo, to 
make reliable predictions of the image separation distribu- 
tion for the subhalo component. 

It is interesting to compute the posterior PDF of K ex t 
to see e xpected environme nts of lensed satellite galaxies, as 
done bv lQguri et alJ (120051) . In FigureEJ we show the PDFs 
for several different image separations. As expected, they are 
strong function of the image separation: Larger-separation 
lenses tend to lie in more dense environments. For instance, 
at 9 < 1" only ~ 10% of lensing by satellite galaxies have 
large external convergence K oxt > 0.1, while at 9 ~ 3" 
more than half of lenses are accompanied by such large ex- 
ternal convergence. Lenses with 9 ~ 6" are dominated by 
those with very large K ox t, which mean that lensing of such 
large image separations occurs only when satellite galaxies 
lie close to the center of their host halos on the projected 
two-dimensional lens plane. 



5.3 Total Distributions 

In this subsection, we put together the results in the previous 
subsections to derive the total image separation of strong 
lenses (eq. 1251 ). 

Figure |H] plots the total image separation distribution 
as well as distributions for the halo and subhalo popula- 
tions. We find that the halo population dominates the lens- 
ing probability at all image separations. However, the con- 
tribution of the subhalo population is also significant: The 
fraction takes the maximum value of ~ 0.2 at around 9 — 1", 
and decreases rapidly at 9 > 3". Therefore, we conclude that 
the subhalo population should not be ignored for the accu- 
rate prediction of the image separation distribution. Since 
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Figure 8. Upper: The total image separation distribution dPt/dd 
(eq. 1251 ) is plotted by a solid line. The contributions of halo and 
subhalo components are also shown by dashed and dash-dotted 
lines, respectively. Lower: The ratio of the distribution for the 
subhalo population to the total distribution, (dP s /d9) / (dP t /dff), 
is plotted as a function of the image separation. 



most strong lenses have image separations of ~ l" — 2", our 
model predicts that 10% — 20% of lenses are produced by 
satellite galaxies and 80% — 90% is caused by central galax- 
ies. 

An advantage of our halo approach is the ability to ob- 
tain further insight into the image separation distribution 
by computing the contribution of the image separation dis- 
tribution from different mass intervals. For this purpose, we 
consider the following three types of halos, "galaxy" (de- 
fined by (host-)halos with masses of M < 10 13 ft _1 M Q ), 
"group" (1O 13 /i _1 M < M < 10 14 A _1 M Q ) , and "clus- 
ter" (1O 14 /i -1 M < M), and explore the contribution of 
each type. We show the result in Figure |5] We find that the 
distributions of group and cluster have tails toward small- 
separations because of subhalo populations. Each type has 
quite similar distributions of the subhalo population, while 
those of the halo population are very different with each 
other. We also integrate the distributions over the image 
separation to estimate the fraction of each component to 
the total lensing probability, which is summarized in Table 
Again, the fractions are very different for halo popula- 
tions, but the differences are much smaller for subhalo pop- 
ulations. As a result, nearly half of lenses which lie in groups 
and clusters are caused by satellite galaxies (subhalos), in 
contrast with lenses in galactic halos which are dominated by 
the halo population. Finally, we compute fractions of lenses 
in lie in groups and clusters and find that they are ~ 14% 
and ~ 4%, respecti vely. These are roughly consistent with 
iKeeton et ail (l200Ct) who predicted that ~ 20% and ~ 3% 
of lenses lie in groups and clusters, respectively (a part of 



Table 1. The fraction of lensing probabilities, integrated over the 
image separation. The entry "total" indicates a sum of lensing 
probabilities for halo and subhalo populations. 



type 


galaxy 


group 


cluster 


([h-^M e ]) 


(< 10 13 ) 


(10 13 - 10 14 ) 


(> 10 14 ) 


halo 


0.74 


0.08 


0.02 


subhalo 


0.08 


0.06 


0.02 


total 


0.82 


0.14 


0.04 



the differences may be ascribed to the different definitions 
of "group" and "cluster"). 

How well does the total image separation distribution 
derived here account for observed image separation distribu- 
tions? We compare our model prediction with the following 
two observed distributions: 

• All 75 lensed quasars discovered to date. This sam- 
ple includes both optical and radio lenses and are quite 
heterogeneous, but it has an advantage of the large num- 
ber of lenses including several intermediate-separation lenses 
(3" < 6 < 7") as well as one large-separation lens whose sep- 
aration is 6 ~ 15". 

• 2 2 radio lenses discovered in the CLASS llMvers et alJ 
2003). Although this sample contains smaller number of 
lenses (and all lenses have 6 < 5"), it is much more homo- 
geneous and has well-defined selection functions. Although 
among the 22 lense s 13 were selected fo r a statistically well- 
defined subsample ijBrowne et al.ll2003u . we use all 22 lenses 
simply to increase statistics. 

The results are shown in Figure 1101 We find that our 
model explains both observed distributions reasonably well. 
An exception is that for all lenses the observed number of 
sub-arcsecond (9 < 1") lenses appears to be smaller than 
expected. However, this is clearly due to a selection effect 
that sub-arcsecond lenses are difficult to find particularly in 
optical lens surveys. Indeed, for the CLASS lenses, which is 
a lens survey in a radio band, the discrepancy is less signifi- 
cant. Thus we conclude that our model well accounts for the 
observed image separation distributions. 

It should be noted that the theoretical predictions is 
based on many assumptions. The most important ingredi- 
ent of our model is the mass-to-luminosity relation defined 
in equations @ and I1UI . To see how the result is dependent 
on the adopted relation, in Figure 1111 we show the depen- 
dence of the total image separation distribution on the mass- 
to-luminosity relations. As seen in the Figure, the shape 
(rather than the overall normalization) of the image separa- 
tion distribution is quite sensitive to the mass-to-luminosity 
relation. This indicates that precise measurements of the 
image separation distribution offer quite useful probe of the 
mass-to-luminosity relation. 



6 SUMMARY AND DISCUSSIONS 

In this paper, we have constructed a model of the image sep- 
aration distribution. Our model is based on dark halos and 
subhalos which are linked to central and satellite galaxies re- 
spectively via simple mass-to-luminosity relation. For dark 
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Figure 9. The contributions of different types of halos on the image separation distribution. We consider the following three types: 
"galaxy" which is defined by (host-)halos with masses of M < 10 13 h~ 1 M Q , "group" by lO 13 h _1 M0 < M < 1O 14 /i _1 M , and "cluster" 
by 10 14 /i — 1 Mq < M. Left: The total (halo plus subhalo populations) image separation for each type is plotted. The sum of distributions 
of three types is shown by thick solid line. Right: For each type, contributions of halo (dashed) and subhalo (dash-dotted) populations 
are shown separately. 
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Figure 10. Observed image separation distributions (histograms) are compared with the theoretical prediction presented in Figurc|S] We 
consider two observed distribut ion: (i) The distribution of all 75 lensed quasars discovered to date, and (ii) that of 22 lenses discovered in 
the CLASS iMvers et alj|2003ft . The observed distributions were shifted vertically to match their normalizations to the theoretical one. 



halos and central galaxies, we have considered baryon cool- 
ing in a dark ma tter halo using an im proved adiabatic con- 
traction model of lGnedin et alJi2004) . For satellite galaxies, 
we have considered the mass associate with their host halo. 
Our primary interest is to quantify the contribution of the 
subhalo population which has been ignored in all previous 
studies. Our results are summarized as follows: 

• We predict that 80% — 90% of lenses should be caused 
by central galaxies (halos) and 10% — 20% of lenses should 



be produced by satellite galaxies (subhalos). The fraction of 
the subhalo population takes the maximum at 9 ~ 1", and 
becomes smaller with increasing image separations. 

• The mass distributions of lensing halos at small (9 < 
3") and large (9 > 10") image separation regions are close 
to those of SIS and NFW, respectively. At 3" < 9 < 10" the 
mass distributions are complicated and cannot be described 
by either SIS or NFW. 

• For both halo and subhalo populations, already at 
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Figure 11. The dependence of the total image separation distribution on the mass-to-luminosity relations L(M) (eq. and L S (M S ) 
(eq. ESI ). Left: We change the overall normalizations of L(M) and L S (M S ) by multiplying numerical factors to them. From upper to 
lower lines, we multiply 10 ' 5 , 10 ' 25 , f, 1CT ' 25 , and 1CT ' 5 . Right: We shift the characteristic mass scale M\\ and m\\ in the mass-to- 
luminosity relations by multiplying the same factors to both Mn and mil (upper lines correspond to the cases that smaller factors are 
multiplied). 



9 ~ 3" the effect of dark halos becomes significant, mak- 
ing a simple SIS approximation worse. This means applying 
isothermal profiles to these lenses results in biased estimates 
of parameters such as velocity dispersions and the Hubble 
constant from time delays. 

• Our model predicts that lensing halos of small sepa- 
ration lenses are heterogeneous: Halo-by-halo differences of 
the shape of rotation curves and the fraction of dark matter 
at image positions are large. This is a natural consequence 
of the scatters of the concentration parameter and the mass- 
to-luminosity relation around their medians. 

• In computing the image separation distribution for the 
subhalo component, it is important to take account of the 
external mass which comes from their host halos. This also 
implies very strong dependence of environments of lensing 
satellite galaxies on image separation distributions. 

• Halo and subhalo populations have quite different con- 
tributions of image separation distributions from different 
(host-)halo mass intervals. We predict that ~ 14% of lenses 
lie in groups, and ~ 4% in clusters. Almost half of lenses in 
groups and clusters are produced by satellite galaxies (sub- 
halos), rather than central galaxies (halos). 

• We have compared our model predictions with observed 
image separation distributions and found that they are in 
reasonable agreements with each other. Since the shape of 
the image separation distribution is rather sensitive to the 
mass-to-luminosity relation, precise measurements of the im- 
age separation distribution in well-defined statistical lens 
samples offer powerful probe of the connection between dark 
halos and galaxies. 

In summary, we have constructed a realistic model 
which predicts lensing probabilities from small to large im- 
age separations. The model is quite useful in understand- 
ing lens populations as a function of the image separation. 



The fraction of lensing by the subhalo population, which has 
been derived for the first time in this paper, is important for 
not only an accurate prediction of the image separation dis- 
tribution but also interpreting results of mass modeling of 
individual lens systems since central and satellite galaxies 
may show rather different lensing characteristics. 

We note that there are several ways to improve our 
modeling. First, we have considered early-type galaxies only. 
Although this can be justified because most of lensing galax- 
ies are early-type, the very existence of lensing by late-type 
galaxies indicates that the correct model needs to take both 
the two galaxy types into account. Since the lensing by 
spiral galaxies has qui te small image separation 6 < 1" 
and also is inefficient llKeeton fc KochanekllT993l . we ex- 
pect that the inclusion of late-type galaxies decreases the 
number of sub-arcsecond lenses caused by central galax- 
ies. In addition, we have neglected the scatters of the re- 
lations between galaxy luminosities, velocity dispersions, 
and effective radii. The scatters could affect the quanti- 
tative results, as the scatter in the mass-to-luminosity re- 
lation is important. The redshift evolutions of relations 
are also neglected. However, if we assume the evolution of 
early type galaxies at z < 1 is purely passive, it implies 
that the mass distribution of lensing halos does not change 
across the redshift. Therefore we expect the effect of red- 
shift evolution is rather small. More importantly, we have as- 
sumed si mple spherical halos . In reality, dark halos are quite 
triaxial Jjing fc Sutol [2002) and the triaxiality enhances 
lensing probabilitie s by a few factors at l a rge-separations 
jQguri et alJ 1200,?! : lOeuri fc Keetonl 120041: iHennawi et"afl 
2005a b). On the other hand, small-separation lensing prob- 
abil ities are hardly affected by the ellipticity of lens galax- 
ies: IHuterer et al.l (|2005) showed that the change is less 
than a few percents unless quasars are very bright. Large- 
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separation lensing probabilities are also sensitive to the in- 
ner slope of dark matter halos (e.g..lKeeton fc MadavJl200li 
IWvithe et aUE)mllTa,ka.ha,shi fc Chihall200ll:lLi Ostrikeil 
2002; Oguril l2002l) . We not, however, that the mass distri- 
bution after baryon cooling is rather insensitive to the inner 
slope of dark halos, thus the fraction of lensing by satellite 
galaxies does not depend on the inner slope very much. Re- 
fining our model by incorporating these is beyond the scope 
of this paper, but is of interest for future studies. 
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